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Amendments to the Drawings: 

The attached sheets of drawings include changes to Figures 2 and 7A-7C. These sheets 
replace the original sheets that included Figures 2 and 7A-7C. In Figure 2, reference numbers 
14, 16, 1 8 and 20 have been deleted, and lead lines have been added for Arm 1 and Arm 2. In 
Figures 7A-7C, reference numbers 14 and 16 have been replaced by reference numbers 14a and 
16a, and reference numbers 14b ad 16b have been added. 
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Annotated Sheet Showing Changes 



REMARKS 



The claims pending in the subject apphcation, claims 1-26, stand rejected under 35 
U.S.C. 1 12. According to the examiner, there is a conflict in what is shown in Figures 1 and 2. 

As discussed below, it is respectfully submitted that Figures 1 and 2 are consistent. 
However, to clarify what is shown in these figures. Figure 2 and Figures 7A-7C have been 
amended. 

Figure 1 of the subject application shows a cantilever beam 14 moxmted parallel to the y- 
axis and fixed to a support 18. The beam 14 bends with negligible extension, and the motion of 
the tip of the beam is substantially in the x direction. Thus, the beam 14 is called the x-arm (Arm 
1). Also in Figure 1, a cantilever beam 16 is mounted parallel to the x-axis and fixed to a support 
20. The beam 16 bends with negUgible extension, and the motion of the tip of this beam is 
substantially in the y direction. Thus, the beam 16 is called the y-arm (Arm 2). 

To simplify the analysis of cantilever beam stiffiiess, it is common to schematically 
represent a cantilever by a linear spring at the end of the beam, oriented perpendicular to the 
cantilever beam, and parallel to the beam bending direction. It is understood that this is an 
approximation which neglects the extension of the beam. This representation is common in 
undergraduate mechanical engineering curricula such as at Figure 1 of Glauser et al.. 
Mechanical and Aeronautical Engineering Senior Laboratory, pp. 1-10. See, the section entitled 
"Modeling a Cantilever Beam." Reference can also be made to page 8 of ASEN 2003: 
Introduction to Dynamics and Systems Spring 2002, pp. 1-19. See, the section entitled "B.3 
Modeling of Vibrations of Cantilever Beams by Elementary Beam Theory." Copies of both of 
these publications are attached. 

Figure 2 schematically represents the equivalent stiffhess of the beam 14, the x-arm, and 
the beam 16, the y-arm. The bending stiffiiess of the beam 14 is substantially along the x 
direction and is represented by a virtual spring element 14a. The high stiffhess which prevents 
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extension of the beam 14 in the y direction is represented by a virtual constraint element 14b. 
The bending stiffness of the beam 16 is substantially along the y direction is represented by a 
virtual spring element 16a. The high stiffhess which prevents extension of the beam 16 in the x 
direction is represented by a virtual constraint element 16b. 

By way of further explanation, a simplified plan view of the system of Figure 1 of the 
subject specification is shown below: 



Top View of Figure 1 




«Arm2 



X Arm 

«Arm 1 ^ 

///>/// ■ 



Combining this view of Figure 1 with that of Figure 2, as illustrated below, shows the schematic 
representation of the equivalent beam stiffhess: 



Figure 2 showino schematic 
representation of equivalent beam 
stirlness. 




/////// 
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Thus, as can be seen from the subject specification, as well as from the above 
explanation, what is depicted in Figures 1 and 2 is consistent with the operation of the claimed 
apparatus for manipulating an object. 

In view of the foregoing, it is submitted that all the claims are now in condition for 
allowance. Accordingly, allowance of the claims at the earliest possible date is requested. 

If prosecution of this application can be assisted by telephone, the Examiner is requested 
to call Applicants' undersigned attorney at (510) 495-3206. 

Please apply any other charges or credits to Deposit Accoimt No. 500388. 




RespectfiiUy submitted, 
BEYER WEAVER & THOMAS, LLP 



Reg. No. 28,411 



P.O. Box 778 

Berkeley, CA 94704-0778 
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Introduction to Experiment 1: Vibrations 

From jackhammers to buildings, automotive 
suspension systems to aircraft gas turbine engines, it 
is obvious to even the most casual observer that 
vibrations are extremely important. In this 
experiment you will examine the vibrations of the 
spring mass, damper system. In part 1 you will look 
at the effect of sinusoidal excitation on a cantilever 
system with small damping. The trust will be to gain 
an understanding of the natural frequency of a 
system. You will vary the mass of the beam, hence 
effecting its frequency response characteristics. In 
part 2 of this lab you will study the effect of 
damping on the frequency response characteristics. 
This part of the lab provides insight into how the 
engineer and scientist can use damping to minimize 
the negative effects of vibrations. Before you begin 
this lab you should review second order differential 
equations with constant coefficients. Keep in mind 
what quantities in the governing second order 
differential equation for the spring mass damper 
system that you are changing throughout the 
experiment. This should be discussed in some detail 
in the lab write-up. 



Part 1: Small Damping 

Required Equipment: 
• Shaker Table 
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• Accelerometer 

• Steel bar 

• Charge Amplifier (I.C.P) 

• Instrumented Hammer 

• Dell Pentium III 800 with National 
Instruments Lab VIEW 5.1 

• Measuring tape or ruler 

• Oscilloscope 

Experiment Apparatus: 

The apparatus consist of a steel cantilever beam 
mounted on a shaker table, which is capable of 
producing frequencies of up to 60 Hz. An 
accelerometer is mounted on the free end of the 
beam, and is used in conjunction with a charge 
amplifier and oscilloscope to monitor the beam's 
response. In addition, extra masses are available that 
can be added to the cantilever to alter its response. 

In this experiment, the oscilloscope is used to 
monitor the electrical signal, which is produced by 
the accelerometer/charge amplifier. 

The oscilloscope trace is actually a graph of voltage 
vs. time, which can be translated or scaled using the 
oscilloscope adjustments. The trace can be 
"stretched" or "shrunk" in either the horizontal or 
vertical directions by adjusting the voltage (vertical) 
or time (horizontal) scales. On the oscilloscope, 
these scales are expressed as either voltage or time 
per division. For example, if the time scale is set to 
10 ms per division, the horizontal dimension of each 
square block on the display grid corresponds to 10 
milliseconds. 

The LabVIEW software will be used to acquire the 
experimental data. The format of this file is found 
following the "Topics for Discussion". 

Elementary Theory: 

A Simple Spring-Mass System 

Many oscillating systems can be modeled as a 
spring-mass system using as differential equation of 
motion. The displacement, y(t), of such systems can 
be found using 
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where m is the mass of the object, c is the damping 
coefficient, k is the spring constant, and F(t) is some 
forcing function. Each term in this expression is 
actually a time-dependent force: my„ is the inertial 
force, cy, is the frictional force, and ky is the spring 
force. Modeling the damping of a system in this way 
assumes that the damping force is proportional to the 
velocity of the mass; this is called viscous damping. 
For a fi-eely vibrating system (F(t)=0), the damping 
coefficient of the system can be found 
experimentally using 



where Yj and Y2 are the values of any two 

consecutive maximum displacements that are one 
cycle apart. 

The natural frequency of a system is the frequency at 
which an undamped system will freely vibrate, and 
can be calculated by using 



Modeling a Cantilever Beam 

The cantilever beam is an example of a system, 
which can be modeled as a simple spring-mass 
system. In order to model the vibration of the 
cantilever beam, the end of the beam is chosen as a 
reference point at which the characteristics and 
response of the beam are measured (see Figure 1). 
An equivalent system is then built that has a 
response, y(t), that is identical to that of the actual 
system. 

The spring constant, k, of the equivalent system is 
identical to that of the cantilever beam, and can be 
calculated quite easily using beam deflection 
formulae. Calculation of the equivalent mass, 
however, is much more difficult because all points 
along the beam,s length do not have the same 
response as the end of the beam. This means that the 




2k 
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equivalent mass, m, can not be determined simply 
by adding the masses nx^^g^j^ and m^^^^, but must be 

found by equating the energies of the two systems as 
they vibrate. This type of analysis is called lumping. 
A discussion of this technique can be found in J. E. 
Shigley, Mechanical Engineering Design^ McGraw 
Hill, 1977, 

y(t) 



M^end 



M^beam 





y(t) ^ Ycos(a>t) 



Figure 1: Modeling the cantilever as a simple spring / mass 

system 

The energy, U, stored in the simple spring-mass 
system is 



2 



where co, is the radial frequency and Y is the Mass,s 
amplitude. To find the energy stored in the 
cantilever, the energies of both the beam and the 
mass at the end of the beam must be considered. 



3 jf-r 



Perfonning the integration, the energy stored in the 
actual cantilever system is found to be 
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140 



Hence, by comparing equations (1) and (2), we see 
that the real cantilever system can be modeled as a 
simple spring-mass system with a mass that is given 



Procedure: 

Measure the dimensions of the beam for use in 
calculation of the theoretical natural frequency. 
Also, record the mass of the extra weights that are to 
be attached to the end of the beam during the 
experiments. The mass of the accelerometer is 25g. 
Be sure to measure only the portion of the beam 
which will be vibrationg. Do not include the portion 
of the beam clamped between the jaws of the vice. 

The Shaker Table Test 

In this experiment, the bars forced frequency is 
investigated. When the forcing frequency matches 
the beam's natural frequency, resonance is observed. 

• Make sure that the bar is clamped securely in 

the vice 

• Turn on the shaker table (Note: the switch for 
the motor is on a small switch box facing the 
wall) 

• Increase the table frequency from 0 to 60 Hz 
without stopping and then proceed back to 0 



• Record the resonance frequency. Allow the 
table to operated at this frequency for only a 
very short time! 

• Turn off the shaker table 

• Repeat the above procedure with each of the 
extra masses attached to the free end of the 
beam. Make sure that each mass is secured 
firmly. 

The Accelerometer Test 

In this test, the natural frequency and damping of the 
bar,s free response is investigated. 



by 



33 




Hz. 
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• Ensure that the LC.P. charge amplifier is 
connected to the accelerometer and to the 
input for channel 1 on the oscilloscope via a 
"T" connector. A separate BNC cable 
connects the interface box to the "T" 
connector. 

• Turn on the I. CP., the oscilloscope and the 
computer. (Note: the amolifier should be 
turned on prior to your arrival by the TA, but 
if not, allow 5-10 minutes for the electronics 
to warm up) 

• Make sure that the needle on the LC.P. unit is 
near the center of the green band. 

• For monitoring purposes, set the oscilloscope 
for the following: 

o Non-storage 
o Internal source 
o Auto trigger off 
o Channel 1 
o A.C. 
o Slope (+) 

• Adjust the vertical position so that the trace is 
on the centerline. 

• Adjust the time and voltage scales so that a 
clean response results when the bar is hit with 
the rubber hammer. Make sure that the end of 
the time scale knob is turned to its "XI" 
setting. 

• Run Lab VIEW VI and follow the instructions 
on the screen. 

Exercises 

• Calculate the theoretical natural frequency of 
the beam with and without the extra masses 
attached. 

• Calculate the damping coefficients for each of 
the cantilever systems. 

• Calculate the natural frequency using the data 
from the accelerometer test. 

• On the same graph, plot the theoretical natural 
frequency vs. mass curve for the cantilever 
beam system and all of the experimental data 
points. 

Topics for Discussion 

• Describe the system,s response (i.e. the 
amplitude and frequency) to different forcing 
frequencies. 

• Compare the various measures of natural 
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frequency. 

• How does the natural frequency vary with 

mass? 

• Discuss the type of damping that is present in 
the cantilever. 



The LabVIEW Vi will generate a tab delimited 
ASCii text document in the following format: 

Vibration Experiment - Cantilevered Beam 
Friday, February 05, 1999 
3:03 PM 

Beam Dimensions: 
Length of Beam (mm): nnnE+n 
Width of Beam (mm): nnnE+n 
Thickness of Beam (mm): nmiE+n 
End Weight (kg): nnnE i n 



System Response (by iteration): 

Frequency (Hz) : Pha se( f) : Vol ta ge( f) : 
nnnEin nnnEin nJinE+n 
nnnE+n nnnE+n nnnE+n 
nnnE+n imnE+n nnnE+n 



nnnE+n 



nnnE+n 



nnnE+n 



Time History (from impact): 
Time (sec): 
nnnE+n 
nnnE+n 
nnnE+n 

nnnE+n 



Voltage (t): 
nnnE+n 

nnnE I n 
nnnE+n 

nnnE+n 



System Response (from impact): 

Frequency (Hz): Phase(f): Voltage(0: 
nnnE t n nnn E i n nnn E f n 

nnnE+n nnnE+n nnnE+n 

nnnE+n nnnE+n nnnE+n 



nnnE+n 



nnnE+n 



nnnE+n 



Part 2: Effects of Damping 
Required Equipment: 

• Model 9210 Drive Harmonic Motion Analyzer 

• 3 cylinders 

• mass disk 

• water 

• glycerin 
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Procedure: 

Before all test, make sure that the analyzer is leveled 
and zeroed. Leveling is important to keep the rod 
from striking the table or rubbing against the 
housing. To level, adjust the feet on the bottom of 
the analyzer. 

After leveling, reset the zero phase reference by [1] 
rotating the pulley on the back of the unit until the 
scale is horizontal and [2] adjusting the height off 
the mass (using the fine adjustment knob at the top 
of the pulley and the gross adjustment by moving 
the tape) to the point v^here the red light (which 
senses the mass amplitude scale) is barely glowing. 
Reset the zero phase reference before each new test. 

Allow the system to come to a steady state before 
recording any data. The results from each test should 
correspond to two steady state equations: 

Forcing function = sin(2fuF^t) 

Mass displacement = ^ ^ ^ 

where A^ is the forcing amplitude, A^^^ is the mass 
ampUtude, is the forcing frequency, and 9 is the 

phase angle. For all test record the following 
parameters. 

• A^- forcing amplitude, which is determined 

by reading the scale on the back of the 
machine. 

• F^ - forcing frequency 

• A^ - mass amplitude 

• - mass period 

• 9 - phase angle (lag: 9 < 0; lead 9 > 0) 

^m' ^^^^ *® digital readout by 

changing the switch. 

Test 1 : Time Response 

Set up the analyzer with an empty cylinder and no 
extra mass. Choose a forcing amplitude and 
frequency of 1 Hz or below so that the rod does no 
hit the bottom of the cylinder and the canes do not 
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come out of the cylinder. After the system comes to 
a steady state, record all of the parameters Usted 
above. Repeat this procedure first with water in the 
cyhnder, then with glycerin making sure to keep the 
some forcing amplitude and frequency. Note that the 
drive can be turned off without adjusting the 
frequency. Record all parameters for each case. 

Add on mass disk to the damping rod (not on the 
vanes) and reset the zero phase reference. With the 
same forcing amplitude and frequency, repeat for 
air, water, and glycerin, recording all parameters. 
(Note: only use one mass) 

Answer the following questions in reference to the 
recorded parameters. 

1. Describe the effects (if any) on the response of 
the system due to the different substances in 
the damping cylinder. 

2. Describe the effects (if any) on the response of 
the system due to the increase in mass. 

3. On two graphs of displacement vs. time, plot 
the forcing function (which is the same for 
both graphs) and all 3 responses for a 
particular mass, using the given equations and 
the values obtained experimentally (not data 
points). 

Test 2: Frequency Response 

Remove the extra mass and reset the zero phase 
reference. By adjusting the frequency, find and 
record the approximate value for the natural 
frequency of the system in air. Please do not le the 
rod bang on the table more than necessary. 

With the same forcing amplitude as before, record 
values for all parameters for six different forcing 
frequencies: three below the natural frequency and 
three above it. Repeat for both water and glycerin. 
(Note: make sure you're not too close to the natural 
frequency) 

Approximate values for natural 
firequency : 

Air: 

Water: 
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Glycerin 

Plot (Aj^/A^) vs. CO and 9 vs. co for each of the three 

damping conditions (two graphs with three plots 
each). 

From these plots, refine your estimate of the natural 
frequency. This is where A^^ tends toward infinity 

for no damping and where 9 = -90° . 
Revised natural fi^equencies: 
Air: 
Water: 
Glycerin 

Replot your results on two graphs, (A^^/A^) vs. co* 

and 9 vs. CO*, where co* = co/co^^. Both of these 

graphs will have six plots ^ make them as neat and 
as large as possible. 

Describe in terms of the recorded parameters what 
effect the forcing fi-equency has on the response of 
the system. Be sure to include the relevant 
uncertainty analysis for the measured quantities (see 
notes in lab manual). 



This site is maintained by Jeffrey Taylor 

Copyright © 1999 Jeffrey Taylor. All rights reserved. 

Revised: August 3, 1999. 
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Appendix A. Natural Vibration Modes and Mode Shapes 

A structural element can vibrate with different frequencies depending upon the nature of excitations. 
For example, when a cantilever beam is placed on the shaker and subjected to varying frequency 
excitations, one observes that at particular excitation frequencies the beam vibrates violently. These 
violent frequencies are called resonant modes or natural vibration modes. One purpose of vibration 
analysis and vibration experiment is to identify the resonant vibration frequencies so that the 
expected operational conditions would not trigger those modes. 



V. Unif arm Beams 

(Tmrisvt'rse or betidiiig v^ibratioDS) 
The aame general formula lioUl* for all the following cws. 
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Table 1 Modes and Mode Shapes of Beam Vibrations 
(Source: Mechanical Vibrations by J, P. Den Hartog, Dover Pub., 1985, p.432) 
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Fig, 1 Selected Vibration Modes and Mode Shapes of a Boeing 777 Plane 

(Source: Dr. Michael Mohaghegh of Boeing Co.) 



Table 1 illustrates the analytical frequency formulas for various boundary conditions anddifferent 
mode shapes. Note that the natural frequency of a uniform is expressed as 



(On = a^ 




(1) 



where a„ is a coeflBcient given in Table 1, £ is Young's modulus of the beam material, / is the 
second area moment which becomes Irea angle = bh^ /\2 for a beam of rectangular cross section 
with its width b and thickness h, and is the mass per unit length, that is, /x = py4 = pbh for a 
rectangular cross section. 

In order to gain insight into the overall vibration characteristics, selected modes and mode shapes 
of a Boeing 777 Airplane are shown in Fig. 1 . Once sufficient natural vibration modes of the model 
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airplane are identified, it is important to understand in what ways the vibration characteristics of 
substructural components contribute to the overall vibrations. To this end, we disassemble the plane 
and carry out component vibration testing of the wing, the fuselage, and the tail boom separately. In 
this experiment, a cantilever beam is used to aid you to gain insight into how component vibration 
modes and mode shapes would influence the overall vibrations. In addition, Figure 2 illustrates first 
five natural frequencies and mode shapes of a VC-10 plane, ranging from 1.85 to 4.28 Hz. Note 
that assembled wings, fiiselage and tail boom may or may not vibrate with the elemental natural 
modes. In addition, the in-flight vibration modes would not be the same as the on-ground airplane 
modes. 



1st antisymmetric mode 
1.85 Hz 



2nd antisymmetric mode 
2.56 Hz 




3rd antisymmetric mode 
2.92 Hz 



4th antisymmetric mode 
3.96 Hz 



5th antisymmetric mode 
4.28 Hz 



Fig. 2 Vibration Modes and Mode Shapes of a VC-10 Plane in flight condition 
(Source: Vibration by R. E. D. Bishop, Cambridge University Press, 1979) 
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APPENDIX B. INTRODUCTION TO VIBRATION THEORY 
B«l Introduction 

VIBRATIONS is a terminology used to describe a broad range of phenomena, both natural and 
man-made, ranging from the unceasing flapping wings of a honeysuckle, the swaying back and 
forth of trees during storm wind, to beautiful violin sounds and the rattling of the steering wheel 
of a car in motion. Our bodies are an embodiment of vibratory phenomena which Prof R. E. B. 
Bishop • eloquently summarized in his 1962 Christmas Lectures at the Royal Institution, London: 
"After all, our hearts beat, our lungs oscillate, we shiver when we are cold, we sometimes snore, we 
can hear and speak because our eardrums and our larynges vibrate. ... We move by oscillating 
our legs. We can not even say 'vibration' properly without the tip of the tongue oscillating. And 
the matter does not end there - far from it. Even the atoms of which we are constituted vibrate." 

In this course, we will primarily focus on the vibrations of mechanical systems, or engineering 
vibrations, that are man-made devices, machines, and transportation systems such as automobiles, 
airplanes and satellites as every successful engineering design must address vibration problems. 
The first design category is to create desirable vibrating conditions. This includes clocks, the 
speed of constant turbine rotations, a concrete mixer, the spin rate of spin-stabilized sateUites, and 
electric massagers, among others. The second design category is to alleviate vibrations altogether. 
This include the bumpy riding quality of a car or a airplane, wind-induced vibrations of bridges, 
excessive noise from washer and dryer, and fatigue failure of critical machine components due to 
prolonged vibrations. It is this vibration alleviation or vibration control that constitute a major task 
to practicing engineers. 

The vibrations of mechanical systems are caused by sudden or continuous disturbances, such as 
aerodynamic forces on the airplanes, oscillations of internal combustion engines in automobiles, 
motors in a washer. Qualitatively speaking, the energy contained in the disturbances are transmitted 
to the mechanical systems and finds its way to propagate throughout the system in a way analogous 
to the way flood water finds its way through the flooded region. Unlike the flooding water that 
flows usually into lower altitude, however, the energy carried by mechanical vibrations from one 
part of the system may reverse its transmission path and reverse its flow path. This can be seen 
when one stretches the spring hanging vertically from the ceiling and let it go for its free vibrations 
to continue. In other words, oscillations of a mechanical system means that the system energy can 
move back and forth. 

Therefore, the designer must assess whether any component of the mechanical system break during 
the initial violent disturbance stages, subsequent steady-state vibrations, degradations of system 
performance during the post-disturbance period, and potential fatigue failure due to prolonged 
vibrations. These four distinct vibratory response characteristics are referred to transient vibrations, 
harmonic oscillations, free vibrations, and fatigue cycles. Sudden gust winds acting on an airplane 
in flight is of transient nature, and the subsequent post-gust vibrations are of free vibrations. Fatigue 
cycles are primarily due to normal operational loads. For example, a woman's heart would have 
vibrated 2.95 billion times assuming 72 beats per minutes for an average longevity of 78 years. A 
mechanical component in an airplane whose fiindamental frequency of 72 Hertz would have vibrated 
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for about 370 million times when the plane has flown about 87000 hours ( 8 hours everyday for 30 
years). A tragic feature of fatigue failures in both humans and machine components is that they 
lead to catastrophes when fatigue failures occur. 

B.2 Frequency of a Single Mass and Spring System 

Consider a single mass m attached to a spring k that is subjected to an applied force / (0 as shown 
in Fig. 1. Applying Newton's second law, one obtains the equation of motion for this system as 

mx(t)+kx(t) = fit) (2) 

where x{t) is measured from the static equilibrium condition under its own weight. Of particular 
interest is the nature of its free vibration, that is, the vibration of the spring-mass system without 
applied force, with an initial displacement jc(0) and/or an initial velocity condition i(0). The 
differential equation to be solved can be written as 

Single Spring-Mass Model 




/=sinu;/ 
Fig. 1 Single DOF Spring-Mass Model 

m x{t) + k x{t) = O, with ;c(0) = xq, x(0) = xq (3) 
We assume the solution of this equation in the form 

x{t) = A cos pt + B sin pt (4) 
where A p and B are constants to be determined. Substituting 

x(t) = -p^(A cos pt + B sin pt ) (5) 

and (4) into (3), we obtain 

(-p^ m+k) (A cos pt + B sin pt) = 0 (6) 
As A and B are not necessarily zero, one must have the term in the parenthesis to be zero: 

i-p^ m+k)=0 =^ p = ^Jkjm {r ad I sec) (7) 
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Substituting this into (4), we obtain 



x(t) = A cos y/k/m t + B sin yfkjm t 



(8) 



The integration constants A and B can be determined by using the initial conditions (x(0) = 
jcq, x(0) = xo) as follows: 



x(0) = xq = A cosO + 5 sinO 
i(0) = jco = y/k/m {-A sinO + B cosO ) 



A—Xq 



B = xo/y/k/m 



(9) 



Finally, substituting this into (3) we obtain 



x{t) = xq cos pt + xo sin pt, p = 



m 



(10) 



Single Degree of Freedom Response due to Unit Initial Spring Extension 




Time 



Fig, 2 Response of Single DOF Spring-Mass System to Unit Initial Displacement Condition 
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Single Degree of Freedom Response due to Unit Initial Velocity 

0.2 1 1 1 1 1 1 1 1 1 r- 



0.15- 



0.05 



.2 
Q 




-0.05 



-0.1 



-0.15 



-0.2 



Fig. 3 Response of Single DOF Spring-Mass System to Unit Initial Velocity Condition 

Figure 2 shows the response x{t) when the spring is extended unit displacement and let to-and-fro 
motions (vibrations) continue. In Fig. 3 we illustrate the response x(0 vs. time when a unit velocity 
is applied. Note that the unit displacement response and the unit velocity response are 90^ out of 
phase. The frequency of the model system is given by 



p J — (rad/sec) 
Y m 

The frequency p in (rad/sec) unit can be converted into to Hz (Hertz) by: 

P =2;r/ => /= A ^Hz) 
Note that, when the mass is fixed, the fi'equency is proportional to 

fcxVk 



(11) 



(12) 



(13) 



In practice, whenever the spring constant is altered, by using thicker beams or heavier springs, the 
mass is also changed. 

An important fact is that as long as the ratio of the spring to the mass k/m is constant, the system 
frequency remains constant. This provides not only a design criterion for vibrations, but also 
provides the possibility of scale-model testing by constructing a laboratory model according to the 
relation: 



7 



ASEN2003: Introduction to Dynamics and Systems Spring 2002 



K k 

p^oc — (x- (14) 
M m 

where (K, M) and (k,m) refer to the real structure and the laboratory scale-model, respectively, 
for which one usually selects M » m. 

B3 Modeling of Vibrations of Cantilever Beams by Elementary Beam Theory 

Consider a cantilever beam that is subjected to a tip load P. If the tip load is suddenly removed, 
the beam will continue to vibrate, exhibiting its free vibration mode as shown in Fig. 4a. 




Simplified Modeling 
from Beam to Spring 
Mass Model 



a) First Vibration Mode Shape of a Cantilever Beam: 
A cantilever beam can be thought of consisting several 
discrete masses as shown above. Note that the tip mass 
velocity is the greatest while that of the mass closest to 
the fixed end is the slowest. Thus, the contribution of 
each mass to the system kinetic energy varies according 
to the vibration mode shape. 




k = 3EI 



0.242 m. 



b) Equivalent Spring-Mass Model 
of First Mode Cantilever Vibration 



Fig. 4 Modeling Cantilever Beam by an Equivalent Spring-Mass System 

The static deflection of the beam, w{x), along the beam length is given by 

P 



yix) = 



6ET 



(x^ - 31x0, 0<x <L 



(15) 



where y, P, E, Z^^, and L are deflection along the beam span x, the tip load, Young's modulus, the 
area moment of inertia of the beam cross section, and the beam length, respectively. Note that the 
beam deflection along the beam span is of a cubic shape. 

At the tip end, we obtain with x = L the following force vs. deflection relation: 



2EI 

P = kSr k=-jf^, S=y(L) 



(16) 



where S is the tip deflection and k is the corresponding spring constant. 

If the beam vibrates according to this particular characteristic deflection shape or commonly des- 
ignated as mode shape, we can express the dynamics of this vibration characteristic or mode 



mSit)+kS(,t) = P, k = 



13 



(17) 
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It turns out that the mass m in the above equation is not equal to the mass of the total cantilever 
beam mb = pAL. This is because the spring-mass equation (17) is an idealization of a cantilever 
beam as shown in Fig. 4b. The determination of the modal mass can be obtained from the system 
kinetic energy equivalence considerations. Note that the kinetic energy of the spring-mass system 
is given by 

T=^mS^{t) (18) 

In order to obtain the modal mass m, we first express the kinetic energy T in terms of the distributed 
velocity of the beam, y(t). We then express y(x, t) in terms of the tip velocity &(x, t). This is 
derived below in a step-by-step manner: 
• First, express from ( 1 6) 

P 38 



Eh. 

Second, substitute this into (IS) to obtain 

8 



(19) 



y{x,t) = :^{x^ -ZLx^) (20) 



Third, differentiate the above equation to obtain the velocity expression 

8 



7(^,0 = ^(^^-31x2) (21) 



Fourth, obtain the kinetic energy of the beam by 



Jo Jo 



(22) 



Fifth, substituting y from (21) into the above equation, we obtain 

r=ijf fpAdx = \ j\^f {x^ -3Lx^fpAdx 
which, when the indicated integration is carried out, yields 



(23) 



] ^ 33p^Z , ^ s2 _ 1 _ j2 
140 



The last expression in the preceding equation is realized by making the total kinetic energy of the 
beam equal to the kinetic energy of the spring-mass system. 

Therefore, the modal mass m is obtained to be 

m = (-f^) = (-^) = 0.2357m6, mt = pAL (25) 
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• Finally, the modal equation (16) becomes 

(26) 



mm+km = P. m = (^^^), ^=^^^- 



140 

The vibration frequency of this modal equation is thus given by 



13 




(27) 



* Equation (24) gives ax = 3.5675 which has about 1.35% error compared with the value obtained by a more 
advanced theory (Mechanical Vibrations by J. P. Den Hartog, pp.429-433). 

Finally, often in a simplified modeling of complex structures by a beam, a concentrated mass is 
introduced along the beam span. For example, fuel tanks along the wing span may be modeled as 
concentrated masses. Figure 5 represents an example of a concentrated mass acting on the tip of a 
cantilever beam. 




Tipmass attached 



Simplified Modeling 
from Beam to Spring 
Mass Model 




Fig. 5 Modeling Cantilever Beam by an Equivalent Spring-Mass System 

The equation of motion for the fundamental mode of the cantilever beam is obtained by adding the 
concentrated mass to the modal mass derived in (26): 



k - 



(28) 



m S{t) + k Sit) = P, m = O.lAlpAL + M, 
Note that the spring constant is the same as before. 

Fundamental-mode (or first mode) of beam vibrations, with different boundary conditions such as 
simply supported and built-in beams as well as with other concentrated masses, can be similarly 
obtained. 

B.4 What Are Resonance Frequencies? 

Suppose you are driving a car on an isolated test road while varying the speed of your car at will. 
The test road is constructed to simulate realistically a variety of road conditions, including back 
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country roads, highways, and freeways. Chances are that you may experience either all or some the 
following conditions: a sudden rattling of the steering wheel, the seat vibrations, and the vertical 
motions of the car body itself If you experience violent vibrations (hopefully luxurious cars don't 
subject you to such unpleasant vibrations for their money's worth!) from your seat, this vibration 
can be modeled as follows (see Fig. 1 of Section 2): 

mx(t) +k x{t) = f sin (ot (29) 

where m is the mass of you and the seat system, k is the spring constant of the seat (and partially the 
elasticity of your body in seated condition), / is the excitation strength being transmitted from the 
car floor to the seat frame, and co is the excitation carrier frequency (or simply excitation frequency). 

Dividing (29) by the mass m, we obtain 



x(t) + x{t) = - sinco t, P =J— (30) 
m y m 

where p is called the natural frequency of the spring-mass system, which should not be confiised 
with the excitation frequency c*;. 

The solution of (30) can be obtained by assuming the following form: 

x(t) = A sinco /, A = constant (31) 
Substituting this into (30), we find 

f m 

-0? A sino) t^' p^ A since> t = — sino; / i-op-A + p^A ) svacot = 0 (32) 

^ f 



Since sin cot does not vanish in general, we obtain 



Since mp^ = m ■ — = k 
m 

Substituting this into the assumed solution form (31) leads to the following important expression: 

Note that the term { in the above equation represents the static deflection x^, from 

kx„ = f ^s, = j (35) 
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Using this relation, (34) can be written as 

— = nico) sin CO /, n(co) = \ (36) 

Observe that, as | sin a;/ 1 < 1, the ratio of the dynamic to the static deflection ^ is characterized 
by the impedance function H{co). 



I 




Frequency Ratio (^) 

Fig. 6 Effect of Vibration Amplitude as Ratio of Excitation and System Frequencies 

Figure 6 illustrates the maximum vibration amplitude as function of the ratio ^ , that is, the excitation 
frequency (co) vs. the natural frequency (/?). Note from Fig. 6 that, when j = 1, the absolute 
magnitude of the system impedance Hico) attains infinity. As the spring cannot be extended to 
infinity, the mass will have to separate away from the spring. This phenomenon is known as 

resonance. 

The foregoing analysis of resonance phenomenon indicates that a good designer must avoid this 
deleterious happening from the outset. Observe that the system frequency should be either lower, 
say, co/p < 0.5 or higher 1 .25 < co/p in order to avoid the violent resonant vibrations. In practice, 
the excitations can vary, suggesting that the natural frequency should be sufficiently away from 
potential excitation frequency ranges. For example, the excitation frequencies that are transmitted 
to your car seat from the car floor would not normally consist of a single frequency. Characterization 
of potential excitation frequencies for specific design components, therefore, constitute a major 
modeling issue. 

It must be reminded that the foregoing single degree of freedom spring-mass system is either an 
idealization or is concerned with one mode of a structural member. For example, the first three 
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vibration modes of a cantilever beam is shown in Table 1 taken from Mechanical Vibrations by J. 
P. Den Hartog. 



V» Uniform Beams 
(TniHavcrac or bOfidiiig vibrftUons) 
The same general formula bolrlft for all the folic wing cw/rb^ 




where EI is the bending eiiffnesi^ of the eeciion, t is the length of the beam, pi is the 
maas per unit iDtigth - IV /gl^ and a« ie a ni]mtir[(»i] cvrulnutj diili^rt'nt for eiuih rose 
and liatud bclnv 



Cant.jiftvfir or "clamped-freo^' beam 
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22.0 
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20U,(I 
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«1 




jri..i 




aa half a ''olamped-ckjunw^" for 




a 






even o^numbeis 


«a 




104. 
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e» 
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Table 1 Modes and Mode Shapes of Beam Vibrations 
(Source: Mechanical Vibrations by J. P. Den Hartog, Dover Pub., 1985, p.432) 



In order to appreciate further the importance of resonance phenomenon, let us assume that the 
excitation frequencies co can be parameterized as 



EI 



where c is a constant parameter. 
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If c ^ 22, then the excitation may not trigger the first mode and third modes. Nevertheless, it would 
certainly cause the beam to vibrate violently with the second mode. Hence, the designer must 
consider all the possible vibration modes and modify the design, if necessary, to avoid resonance 
catastrophes. 
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B.5 Spring-Mass System with Damper 



A well-designed car seat dampens out a wide range of excitations. In addition, there is friction 
between the seat and the passenger. To model such damping phenomenon, we modify the spring- 
mass system shown in Fig. 1 by attaching a damper c as shown in Fig. 7. 

Single DOF Spring-Mass-Damper Model 




X 



, f. 

Figure 7 Single DOF Spring-Mass-Damper System 

The equation of motion for this system is obtained by modifying that of the undamped system (28) 
to yield: 

m x(t) + cx(t)+kx(t) = f sin cot (38) 
It turns out that the solution of the above damped vibration equation can be obtained by 

x{t) = A sin cot + B cos cot (39) 

It is recalled from (30) that one needed only the first term of (39) for solving the undamped system 
(28). 



Substituting 

x{t) =^ co(A cos cot ~ B sin cot) 
x(t) = o/(A sin cot + B cos cot) 

into (38), one obtains the following equation: 



(40) 



(-co^mA - cocB + kA) sincot + {-co^mB + cocA + kB) cos cot = fsincot (41) 

By equating the lefthand and the right-hand side of the same coefficients of sinct;^ and cos cot, one 
obtains 



(-co^mA - cocB +kA) = / 
i-co^mB +cocA+kB) =0 



(42) 
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Equation (42) can be expressed as a 2-by-2 matrix equation: 



k — co^m —coc 
(oc k — aP'm 



(43) 



In order to simply the resulting solution, we divide the first and second rows of this equation by k 
to obtain: 





' A 




f ^ 1 

k 




B 




0 



Two further substitutions are made to the above matrix. First, we note that 

m 1 1 



(44) 



k 1 p2' 

m ^ 



k p 



Second, we parameterize c by 
so that we have 



c = 2^mp, 0 < f 



c m c m 1 CO 

co~r = -(o-r = 2^pco— = 2f ~ 
m k m k p^ p 



(45) 

(46) 
(47) 



Using (45) and (47), (44) is simplified to 



Now, solving for {A, B) we obtain 

A = 

B = 





A 




\ f 1 

k 




B 




0 




f 





(48) 



[l-(f)2]2 + 4f2(p2 k 



(49) 



^ p 



_ / 



Finally, if we substitute {A,B) into (3 8), the resulting solution becomes algebraically quite complex. 
Therefore, we simplify (38) by using a trigonometric identity as follows. Note that (37) can be 
modified as 



x{t) = A sin a;/ + 5cos(w/ = yj + B'^ • (- 



V^2 + ^2 

= ^ A^ + * (sina;/cos<^ + sin0cosa>O 

y B 

= y/A^ + B^ ' sin{o)t +(t>), tan0 = — 

A 



sin o)t + 



B 



+ 52 



COSO)/) 



(50) 
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Substituting (49) into the above expression, we obtain the solution of the spring-mass-damper 
system (37) under a sinusoidal excitation as 

xit) = , • sin(££;/ + 0), tan 0 = 

y[l_(p2p + 4^2(p2 [1 - (p2] 



(51) 



x{t) _ x(,t) 



= n sin(cot + 0), n = 



y[l-(f)'P + 4^2(^2 




Frequency Ratio (— ) 

Figure 8 Amplification Factor of Single DOF Spring-Mass-Damper System for Harmonic Excitations 

The parameter f is called a critical damping ratio and affects the vibration responses as follows: 

if f z= 0.0, the system is undamped (see Eq, (35) ) and continues to oscillate; 
if 0 < < 1.0, the system is damped but eventually reaches to its static displacement; 
if ? > 1 .0, the system is critically damped and does not oscillate. 

(52) 

For a given damping ratio (O, the amplification factor Ti is uniquely determined by the frequency 
ratio (u;/p). This is illustrated in Fig. 8 for various damping ratios. For the undamped case 
(f = 0.0), the result is the same as the one abeady shown in Fig. 6. 
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Figure 8 indicates that an ideal choice of the seat design parameters (the seat-and-passenger mass 
m and the spring constant ks) is to choose p = to be several times smaller than the largest 
potential excitation frequency c*;. However, this usually implies that the seat spring constant ks to 
be impractically small; a heavier passenger may cause the seat frame to buckle! This is where the 
seat damping can play a major role in mitigating vibration as shown in Fig. 8. Typically, if the seat 
damping characteristic is chosen such that the resulting damping ratio becomes 

0.707 < i: < 1.0 (53) 

it is considered a good compromise (why?). This and additional features to be discussed later in 
Class Notes may serve as a basic design guide for vibration design. 

B,6 Predicting the Vibrations of a Real Structure from Laboratory Experiments 

Suppose you are tasked to design a scaled-down model structure, carry out a laboratory test of its 
vibrations, and predict the vibrations of a real structure. This is often necessary both for cost, rapid 
design cycle and facility considerations. For example, it may be feasible to design several full-scale 
car frames and test them for their vibration characteristics. However, it would not be feasible 
to design and manufacture several full-scale wings or fuselages, and test them for their vibration 
characterizations. Table 2 lists a typical scaling relations for correlating a laboratory-model car 
crash test to a full-scale car crash. 



Table 2: Scaling of Laboratory vs. FuU-Scale Design Variables 
(5>1: Scale Factor) 



Variables 


Full Scale 


Test System 


Length 


L 


e = L/s 


Displacement 


X 


x = X/S 


Time 


T 


t = T/S 


Velocity 


V 


V 


Stress Intensity Factor 


Kc 


kc = Kc/VS 


Stress 


a 


a 


Strain 


€ 


€ 


Density 


P 


P 


Young's Modulus 


E 


E 



For the purpose of correlating typical laboratory vibration test data to a full-scale structure, unfortu- 
nately, there has not been a straightforward relation for us to employ. One reason for this difficulty 
is that the frequency equations for plates, shells and cylinder-like structures take on different forms 
than for the case of beams. By restricting ourselves to the vibrations of beam or beam-like structures, 
we may utilize the following relation: 



18 



ASEN2003: Introduction to Dynamics and Systems Spring 2002 



-4 = TIJT (54) 

where the subscripts / and /refer to the full-scale and laboratory models, respectively. Nevertheless, 
many airplane vibration modes can be characterized in terms of simplified beam models. 
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